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Abstract 

Rules for quantizing the walker+coin parts of a classical random walk are provided by treating 
them as interacting quantum systems. A quantum optical random walk (QORW), is introduced by 
means of a new rule that treats quantum or classical noise affecting the coin's state, as sources of 
quantization. The long term asymptotic statistics of QORW walker's position that shows enhanced 
diffusion rates as compared to classical case, is exactly solved. A quantum optical cavity imple- 
mentation of the walk provides the framework for quantum simulation of its asymptotic statistics. 
The simulation utilizes interacting two- level atoms and/or laser randomly pulsating fields with 
fluctuating parameters. 
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I. INTRODUCTION 



Discrete quantum random walks on a line, (for a review in a more general 

setting see [l|) is a paradigmatic construction of a quantum system performing a motion 
similar to the usual classical discrete random walk on a line, but with a number of important 
differences like the quantum treatment of its coin and walker systems, the role of quantum 
entanglement, novel diffusion and hitting properties of its motion etc. The asymptotic be- 



havior of such walks has been studied 



ied in^. 



10j,[2l|. Also important theoretical suggestions 



about the utilization of those properties in the construction of quantum algorithms that 
would outperform classical rivals in tasks such as searching a database etc, have been put 
forward, see e.g. {4]. The effect of decoherence on the evolution of discrete quantum walks 
has been studied in e.g. [l^,[l6|. Experimental proposals also exists concerning e.g. quan- 
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12|, 13|. The coin control 



tum coin-tossing, and quantum diffusion see e.g. [ll|,[26|, 
over discrete quantum walks on graphs has been studied in [20(]. Continuous time quantum 
walks on graphs, which are defined without the use of a coin, were studied in e.g. , and 

their potential implementation for the construction of fast search algorithms were studied 
in 0. 

The present work extends previous ones on the so called V'' models of quantum discrete 



walks on a line 



22 



231] 2J] 251]. Its aim is twofold: first to argue that a general and systematic 



framework of quantization of classical walks is possible and well physically motivated, and 
second to study the walk beyond transient dynamics in its asymptotic regime and to show 
the possibility of simulating its statistical behavior in terms of another quantum system, 
introducing thus the quantum simulation of a quantum walk. This double aim is achieved 
by introducing the Quantum Optical Random Walk (QORW); a mathematical model that 
refers to a typical cavity QED type of arrangement where a beam of two-level atoms passes 
through an optical cavity. [Technical caveat: although the mathematical description of 
cavity field done here is not the usual one i.e. in terms of the boson degree of freedom, but 
instead is in terms of states with positive and negative valued energy also (canonical algebra 
vs. Euclidean algebra, see below), appropriate truncations can be introduced in order to 
have only positive energy states (c.f. e.g. [12I); this remedy for the problem however will 
not be discussed here. ] 

The outline of the paper goes as follows: in Chapter II, the study of quantization rules is 
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developed, the QORW is introduce in Chapter III and several example of walks are solved; 
Chapter IV treats the walk in its asymptotic regime and introduces the technique of quantum 
simulation, examples of numerically simulated asymptotic probability density functions are 
provided; conclusions are summarized in last Chapter V. 

II. QUANTIZATION RULES FOR CLASSICAL RANDOM WALKS 

The essential feature of a simple QRW on a hne is the promotion of mathematical 
correspondence: left/right^ head/tails, between walker^s move directions and coin^s two 
sides, to a dynamical interaction among two physical systems. This is realized by in- 
troducing state Hilbert spaces = span{\m))m£Z and He = span{\+), |— )), for quan- 
tum walker and coin systems respectively. In operates the Euclidean algebra with 
generators the step operators Ej-\m) = |m ± l),(and their their powers e.g. for a > 0, 
{E±)°'\m) = ^±a|m) — \m± a)), and the position operator L\m) — m\m). They satisfy the 
commutation relations [L, E±\ — ±E±, [E+,E-] — 0. Also important is the Fourier basis 
Hyj — span{\4>) — ^ Ylimez < (j) < 27r), which is the eigenbasis of step operators 

viz. E±\(p) = e*"^|(/)). In the coin space the projection operators P± = |±)(±|, are needed in 
order to realize the coin-tossing that drives the walk. Indeed one step of classical random 
walk (CRW), is described by means of the unitary Vd = P+<^ E+ -|- P_ ® E_. Its conditional 
action on walker states realizes the coin-tossing and the subsequent move of walker. 

Explicitly the coin pc and walker p^, density matrices initially taken in product form 
Pc ® Pw-i are assumed to interact unitarily by the transformation Pc <S> Pw Vci{Pc ® 
p^)T4;. Subsequently the two systems are decoupled by an unconditional measurement 
of the coin subsystem realized by means of the partial trace i.e. Vci{pc <S) pw)Vd 
Trc{Vci{pc <8) Pw)Vci) ; the latter constitues the dynamical realization of the coin tossing 
process. The resulting walker density matrix ey^iiPw) — Trc{Vci{pc 'S> Pu,)K;), written in 
the eigenbasis states of position operator L (see above), provides by means of its diago- 
nal elements the occupation probability distribution p^ — {m\ev^i{pw)\n^) , of the states of 
walker system. This distribution in the course of time steps of the walk n — 1,2,..., is 
identified with the occupation probabilities of the classical random walk on integers i.e. 
Prn = {m\ev^i{pw^)\m) ,m = 0, ±1, ±2, ...,with bias determined by the elements of pc- Due 
to this the Sy^i map is conceived as realization of CRW, which we next seek to quantize. 
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Quantization of CRW is conceived as the incorporation in coin space of an additional 
unitary operator U, the coin reshuffling matrix, so that the one-step operator now be- 
comes V = VciU ® 1. To facihtate conceptual comparisons this procedure was christened 



[/—quantization in ref. 



251 ] ■ Almost all work that has been done in the area of quan- 



tum random walks, has been based on the scheme of quantization or modifications 
thereof. One particular class of such [/—quantized walks are the V^'^— models for which 
the one-step evolution of the walker's density matrix is given by the CPTP map eyk as 
ey,{p^) = Tr,y\p, ® p^)V^^m. 

Next we introduce a generalized version of the previous quantization method, the 
e— quantization rule, which employs a positive and completely positive trace preserving map 
e, acting on the coin density matrices, which is not necessarily taken to be unital, namely 
e{\l) ^ \l. Then the one-step evolution of walker's density matrix of an e— quantized model 
of a classical walk is defined as 



ev{p^) = Tr,V,ie ®l{p,® p^)Vl 



In general for a V'' quantum walk model we will have that 



(1) 



(2) 



where in general a different quantizing CP map e can be used between coin+walker inter- 
actions. To appreciate the changes brought about by 5— quantization we use the adjoint 
action of e.g operator X, on a density matrix defined as AdX{p) = X pX'^ ,with property 
AdXY{p) = AdXAdY{p). We assume that the e map is determined by a set of Kraus op- 
erators as £{pc) = J^i'^iPcSj, or in terms of the adjoint action e{pc) = J^i^^'^iiPc), 2h\- 
To contrast deference between the usual [/—quantization and the proposed e— quantization 
rules we write the evolution map of the former as 



evk{p^) = Tr^{AdVciAdUf {pc® p^) 
= Tr,{Ad{V,iU)fip,(^p^), 



(3) 



and the evolution map of the latter as 



£yfc(Pu;) = AdVciAdSj^ {pc® Pw) 

= Tr, (^Ad{V,iSi)^ {pc®Pu,). 



(4) 



It is evident from above that the e— quantization rule apphes a sum of adjoints at each 
factor of the /c— fold product. This sum in addition to Vd, is determined by the Kraus 
generators of the quantization map e. In the simplest case where there is only one single 
unitary Kraus generator, the e— rule reduces to the [/—quantization rule. The inclusion of 
Kraus generators in the quantization of classical walk, may stem from the fact of a hidden 
quantum interaction between the coin system and another unobserved quantum system, or 
it may be due to some classically fluctuating parametric variance of the coin system, that 
is treated on the average; generically we may say that it is due to quantum or classical 
noise, (see also previous studies of non-unitary models H, 12, 1^, 17, 18, 19|, The number 
of Kraus generators as well as their matrix type provide in any case a flexible framework 
in which the quantization map may incorporate complex physical processes in coin systems 
that may take place within a evolution step of the walk. 



III. QUANTUM OPTICAL RANDOM WALK 

As an application of the e— quantization rule we introduce a new kind of quantum walk. 
To this end we need to introduce a continuous family CPTP maps E = {t ^ e^, t > 0}, acting 
on the space of coin density matrices, where variable t is taken to stand for time. Among 
members of family E, there exists a semi-group composition law, i.e St-^ o et^ = efj+tg-Also 
the identity map Eq = id, is included in E. 

Now we come to the Quantum Optical Random Walk; it physically outlines the crossing 
of an beam of two-level atoms, the coins, through a quantum optical cavity which sustains 
a standing quantum mode identified with the walker system. The walker-|-coin interaction 
realizes a V"^ QRW model, and it also takes into account the interaction of coin with some 
external environment, formalized as some CPTP time dependent map St- This external 
interaction of the coin is taken to have been initiated at some past time 0, and to con- 
tinuously happen in time while atom crosses the cavity. Entering into the cavity at some 
time t, in state £t{Pc), the atom interacts instantaneously, being the coin part of QRW, with 
the walker/cavity mode. For the V"^ model two such co in -|- walker interactions occur: one 
interaction at time t, that changes their combine state as St{pc) ® Pw ^ yd{£t{pc) ® Pw)V^i, 
and a second one later on at time t + that effects the change of state: Vci{£t{Pc) ® Pw)V^i 

Vci{St ® l(K;«(^t(Pc) ® Pw)V^i))V^i. The Vd 's are realized by sudden on-off switches of 
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the mode+atom interaction, possibly by driving the cavity mode off resonance. Then the 
two interactions taken together constitute one step of the walk. Subsequently the atom is 
considered leaving the cavity, the time clock is been reset, and a new atom is entering the 
cavity. (A realistic cavity QED model for implementing QRW in atom+mode interactions 
has been proposed in [12l].) Explicitly the total change of walker's density matrix between 
two successive steps is given by the equation 

pI"^ = £v2(pI""'^) = Tr, [AdV,i ■Sr(E)l- AdV^i ■ ® 1 ■ (pe ® pI""'^)] , (5) 
for n = 1, 2, ...,or more explicitly 

pI"^ = ev^iP^:-'^) = E (^k*(Pc)b)(^k^(K)01)l^)^«+^P^'^^.+/^- (6) 

ijk=±l 

This shows that the e— quantized V"^ walk proceeds with steps of length and 2, on the 
ladder of walker states, with weights determined by the time dependent CPTP e^.This same 
map actually serves as the source of quantization of classical V^i walk. It is important 
however to emphasize that in this walk the physical origin of quantization is not an ad 
hoc imposed unitary rotation matrix in coin space, as it has been in most cases following 
the [/—quantization rule, but instead it is the physical process of natural interaction of 
coin system with some external agent. As an example we can consider the case of a Rabi 
oscillating two-level atom that decays spontaneously. Such physical conditions render the 
£— generalization of quantization for a classical walk a well motivated one. 

The one-step map £y2 previously introduced is determined by parameters; these are some 
e.g A parameter measuring the strength of Et, and the parameters t and r, determining the 
time intervals of coin-|-walker interactions. Those time parameters should be finely tuned, so 
that is possible for the two interactions to take place during the time the atom/coin spends 
in the cavity. This can be decided by selecting the velocity of the atomic beam crossing the 
cavity. 

To proceed with the problem of time evolution we introduce walker's density matrix, 
Pw = Jq^ p{.4>: (j) )\(j)) {(j) \d(j)d(j)' . Due to linearity of evolution we only need to compute 

^ (7) 

then 

£^(P«.) = / V (8) 
JO Jo 
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In formulas above matrix Vd{(f>), is defined in the eigenbasis of step operators as Vd = 
J Vci{(t>)\(t>) {(t>\d(f), where Vci{(/)) =diag(e*'^, e"*"^). Namely in 0-basis we have that E± = e^**, 
where $ = / is a Hermitean phase angle operator. 

The discrete distribution determining the m-th site occupation probability on the walker's 
ladder after n steps is 

1 /'27r /'27r ^ 

Pj^) ^ {m\eUpJ\m) = j-^ 0')e-™(^-<^ ). (9) 

Change to variables 4>± = (j) ± (j) , in the preceding expression shows that if ^4(0+, is 
independent from variable 0+, then the 0+ integral can be carried out and the probability 
becomes 

Pir) ^ {m\el,{p^)\m) = 1- d<f>^ /l"(0_)e-^'"<^- , (10) 

then the can be identified with the characteristic function of the transition probability 

function of a classical random walk[33]. This is a general result independent of quantization 
rule that can be used as a criterion of classicality of quantum walk. Let us give some 
examples of classical walks quantized according to U and e rules, and impose on them the 
criterion of classicality. 

Examples from ?7— quantization rule: 

i) The V model with evolution ev{piv^) = TrcVpc ® Pw ^^V'^, gives Av{(f)+,(f)~.) = 
(+|Pc|+)e*'^~ + (— |pc|~)e~*'^" .This fulfills the classicality criterion and essentially leads to 
classical occupation probabilities drawn from the diagonal elements of walker's density ma- 
trix. 

a) The V'^ model with evolution ev{pw^) = TvcV^ pc® p^'^^V^^ initial pc = |+)(+|, and 
reshuffling matrix [/| = e*5°^^ , been a ^-rotation. This gives I2J] 

y4y2(0+, 0_) = COS^ 0_ — 2 COS 0+ siu 0_. (11) 

The asymptotic characteristic function (cf. eq.( l22l) ). that determines the limiting proba- 
bilities of the walk, is for this example 

/i(0) = -2 [d^A{<P, <p%,^^ = - cos 20. (12) 
Examples from e— quantization rule: 
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in) Assume initially we have a coin being in some mixed state pc =diag(g, 1 — q), 0<q < 
1, and that it suffers spontaneous de-excitation of its upper state with rate A, then its state 
is described by the map [28|, 



etiPc) = So{t)p,Sl{t) + S,{t)p,Sl{t), (13) 
where cos(At) = Vl — e"^-^*, and 

/ cos(At) o\ / o\ 

So{t) ={ , S,{t) =1 . (14) 

y 1 J I sin(At) / 

Such a coin enters the cavity in state £t(pc) =diag(gcos^ At, — gcos^ Xt + 1). Then we obtain 
that 

Av2{<j)+,(j)_) = e-'^'^-{l - qcos"^ Xt) + e'^'^-qcos^ Xtcos^ Xt + qcos^ Xtsin'^ Xt (15) 

For q = |, we have the initial density matrix = |l, as special case. The criterion of 
classicality is also fulfilled here, so the diagonal elements of evolved walker's density matrix 
can be identified with a classical distribution. 

Statement : Direct calculation verifies the general statement that classicality prevails in 
any e— quantized V'' model, independently from the initial condition of coin system, as long 
as the quantization proceeds by employing an e quantization map with Kraus generators 
being matrices having only diagonal or only anti-diagonal elements; such examples of e 
maps are e.g the X, Y, Z, the depolarization channels, the mentioned spontaneous emission 
channel, as well as the transformation induced on the coin/atom after tracing out the bosonic 



field freedom in Jaynes-Cummings model and in its various modifications (c.f. [2J, |25l]). 

iv) If initially we have a coin being in some mixed state pc =diag(g, 1 — q), 0<q < 1, 
but 7^ |, and if the quantization proceeds by using the map 

eipc) = RoPcRI + RipcRl (16) 

where 

^0 = ^1, ^i = ^t/f, (17) 

we then obtain that 



For q = 0,1, we have the initial density matrix pc = | — )(— |, |+)(+|, respectively. By 
comparing examples ii) and iv) we see that we pass from the [/—quantization rule of the 
former to the e— quantization rule of the latter by simply mixing the reshuffling unitary 
matrix U]l with the unit matrix. Such mixing is relevant in cases where there is probabilistic 
uncertainty as to whether reshuffling matrix is applied or not. 

The criterion of classicality is not satisfied here, therefore we have a genuine e-quantized 
random walk. By virtue of eq. ffTSj) . the asymptotic characteristic function is (c.f. eq. fl22l) ) 

M0) = -^ + ^^cos2(/>. (19) 

Closing this section we note that both quantization rules give non trivial models for quantum 
walks and that in some cases, as e.g in the last example the e— rule is presented as a necessary 
amendment of [/-rule when the latter can not be applied since the reshuffling matrix is either 
not exactly known or is not accurately applied. 



IV. QUANTUM SIMULATION OF ASYMPTOTICS 

Asymptotics: The dynamics of quantum walk can be described by the quantum statistical 
moments of the observable of position operator L, evaluated e.g after n steps. We obtain 
for its statistical moments 

27r 

{L% ^ Tr{L^p^:^) = ^Jd<f> [d; [p(0, 0')]]^,=^ • (20) 



Next we study the asjTuptotic behavior of moments when the number of steps n is large. 
In this case we have that 



n 
2^ 





s-1 



,0)M0)^ + OK-^). (21) 







Here we have introduced the asymptotic characteristic function (acf) h{(p), of the walk as 

M0) = -^[5<^A(0,0O],.=.- (22) 
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The scaled by time limiting statistical moment of position operator is then obtained to 

be 



27r 

LY\ 1 



n j 1^ 27r 



#p(0,0)/i(0)^ (23) 







For the {/-quantized Syk model since h[(^^ =ImTrc(y'^^(0)\^'^(0)pc)) as can be obtained 
by elaborating on last equation, we have that 

h{4>) = Tr, [{a + V\(P)aV{(l>) + ■■■ + V^'-\^)aV''\^))p,] , (24) 

where a = U'^a^U is a rotated a-^ Pauli matrix. 

It is important to note that the limiting positional moments provide all necessary infor- 
mation for a sufficient understanding of the ensuing walker asymptotic dynamics, and as 
seen from above these moments are expressed as classical mean values of the powers of func- 
tion h{(f)) of the stochastic variable 0, that takes values around a circle with pdf ■^p{,4', 0)- 
Hence we will seek the exact knowledge of function h{(j)) next, by the technique of quantum 
simulating the system of walk. 

Quantum Simulation: The concept of quantum simulation of dynamics or statistics of a 
quantum system by another quantum system constitute a final goal for Quantum Information 
Science, since an alleged universal quantum computer device would function as an efficient 
simulator of any quantum process. In more modest claims a special purpose quantum system 
could be constructed and set up to interact with its environment, so that its dynamical or 
statistical behavior would simulate the respective dynamics or statistics of a given quantum 
system. The simulator system is possibly different from the original system, in e.g its 
dimension, type of interactions or necessary physical and computational resources required 
from its time evolution. However both original and simulator systems are both governed by 



laws of quantum mechanics |29l. I30!,l3l|,l32|. 



To construct a quantum simulation of asymptotic behavior of [/-quantized walk, we 
proceed by simulating quantum mechanically its asymptotic characteristic function h{(j)). 
Let us refer to eq.( !24l) . and introduce firstly the CPTP map 



= + yK^)<yV{^) + ■■■ + V^'-\ct^)aV'-\ct^)), (25) 

then we cast acf of eq. (l2^ . in the form /i(0) = kTrc{e*^{a)pc)- Equivalently we can express 
acf in the form /i(0) = kTrc{e^{pc)(j), where utilizing the cyclic property of trace, the dual 
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map of the preceding map 5^, has been used that reads 

e^{Pc) = l{pc + V{ci>)p^^\<p) + ■ ■ ■ + V^~\cP)p^^^''-\c|>)) . (26) 

It is now possible to express the scaled s-th moment of quantum variable L, after the n-th 
step of the walk, in the suggestive form 

" 

Let us first elaborate on the first moment taken for s = 1 in last equation; to this end 
we introduce the 0— average of the transformed density matrix 6^{pc), with respect to the 
probability distribution function (pdf), ( ^^2^^ , < < 27r), that reads 

-e(p^) = j'-^eM)- (28) 



Then we obtain the first moment for n » 1, in the form 

This is interpreted as saying that (c.f. first eq. above) in the asymptotic regime of the 
walk, the expectation value of scaled variable L is proportional to the expectation value of 
observable a, evaluated with the coin density matrix been transformed by the 0— average of 
e^p. In the second equation above, it is emphasized that moment {{^))^ is initially defined 
as expectation value of walker's space observable in the limit of large number of steps. 

Next we continue our elaboration with the case for higher moments i.e s > 1. Referring 
to eq. (127ll . the integrand (Trc[o'e(j,{pc)]Y , by means of the property of trace Tr{A ^ B) = 
TrA- TrB, or TrA'^^ = (TrA)"-, is expressed as (Trc[ae^{pc)]y = Trc[cre^{pc) ^ae^{pc)-.- ® 
cr£(^(pc)]- Further use of the property A ® B ■ C ® D = AC ® BD, (dot sign emphasizes 
ordinary matrix product), and of notation A*^" = A^ A^ ... ^ A, for n-fold tensor product, 
yields that 

{Tr,[ae^{p,)]y = Tr,[a^' ■ e^{p,) ® e^{p,) ® ... ® £^(p,)]. (30) 

We also need to introduce the 0— average of the products of asymptotic density matrix 
£(t>{Pc)i appearing above i.e 

^(pf ) = / ^^^^MPc) ® ^0(Pc) ® ... ® e^{p,)]. (31) 
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Combining eqs. fl30|31l) . we cast eq. (l27|) of the asymptotic s-th moment in following form 



(32) 



Preceding equation (!32l) and its s = 1 version in eq. (!29l) . provide a framework for a quantum 
simulation of asymptotic statistics of quantum walk on a line. Such a framework identifies 
the walker system of a QRW with the system that is simulated by a second quantum system 
the simulator which in this case can be identified with the coin system of the QRW in 
question. 

Indeed let us assume that we set up some appropriate Hamiltonian dynamics for a com- 
posite system comprised by a two-level atom identified with the coin system, and an ancil- 
lary system, so that after decoupling the two, by tracing out the ancilla system, the coins 
finds themselves in state B^{pc)- Then the quantum mean value of the observable cr®'', i.e. 
(cr®'') = Tra^''^e^{pc), will be equal with the asymptotic s-th moment, of the walker sys- 
tem of the simulated QRW, by virtue of eq. (l32|) . The problem of constructing appropriate 

Hamiltonian dynamics has no unique solution, so next we provide a solution for the simplest 
nontrivial case of = 2. For this case eq.( l26l) . becomes 

eM = l{Pc + Vi<l))PcV\<j>)). (33) 

To unitarize this transformation of coin density matrix we introduce a 2D auxiliary system, 
where pa = |+)(+| is taken to be its density matrix. Then we find that 

e^Pc) = TraWmp, ® Pc)W{(Py, (34) 

where the unitary matrix W{<p), is chosen to be 

1 , 1 [ 1 V((h)^\ 

W{<P) = ^{l0l+\-){+\®V{<P)-\+){-\0Vm = -=l . (35) 

V2 V2 y-V{(j)) 1 J 

Writing W{<p) = exp H{(j)), the associated Hamiltonian matrix reads H{(f)) = ^[(t+ C?) V(</>)^ — 
cr_®V^(0)], so that if we choose the reshuffling matrix to be a ^-rotation matrix, i.e. V{(f)) = 
Vci{4>)U = e^'^'^^UiL = e*'^'^3g«f 0-2^ then the Hamiltonian suggests a coupling of two spins and 
becomes 

H{(j)) = j[a+0 e-'^'^'e-''*"'' - a ® e^'^'^'e'^"']. (36) 
12 



Having constructed the Hamiltonian interaction that yields density matrix e^{pc), the 
0— averaged matrix e^{pc), written in the form 

^{Pc) = Tr^ I ^^M^ [w{<pr iPa ® Per w{<pr^] . (3?) 



requires the s-fold product unitary W{(l))®^ = e^^ , with associated Hamiltonian the follow- 
ing multiple spin matrix AH = HCi?)l®...®l + ... + l®...l^H. Explicitly it involves 2s 
two-level systems coupled in nearest neighbor form, 

AH = V{(f)y - a_ ® V{(f))] ® 1 ® ... ® 1 + ^1 ® [a+ ® V{(py - a. ® V{(p)] ® 1 ® ... ® 

+^i(g)i® ...o [(T+®y((/))^ - (7.0 1/(0)]. (3 

For example if initially pw = |0)(0|, i.e. p{(p,4>) = 1, then the density matrix of eq. fl371) is 
obtained by a uniform average of unitary similarity transformation generated by the Hamil- 
tonian of last equation. Having constructed e'^^pc), by means of the evolution of simulator 
system, the moments of quantum walker's system are obtained by averaging measurements 
of observable o"®'^, in accordance with eq. (l32!) . 

Classical Stochastic Simulation : An alternative way to implement physically the transfor- 
mation outlined in eq. fl28|) . is to consider the ensemble average of an appropriate stochastic 
unitary rotation acting on a two-level system. In more concrete terms, let us consider the 
angular random variable (rv) 0, distributed on circle by the pdf (p ~ ( ^^2^^ ' < < 2n), 
and the independent discrete rv ly, uniformly distributed over the first k natural numbers i.e. 
1/ ~ ({0, 1, k — 1}). Then we form a transformation of the coin density matrix as follows 
Pc V {(py pcV {(pY'^ . This is a random similarity transformation of the density matrix with 
the random unitary matrix V{(j)y] its randomness is both due to angle (p and due to expo- 
nent ly. Taking the statistical average of this transformation over its two rv's with respect 
to their corresponding pdf's to be denoted by {V {(/)y pcV {(py'^') ^ ^ , we write that 

(y(0)V.y(0)t^>^_^ = I Jpi^^ (j2 V{xrPcV{xyA ^ e{p,). (39) 

This is identical with eq. (!28|) : note also that due to the statistical independence of and 
u variables, the double statistical mean of rotations is obtained by evaluating successively 
the mean for each variable i.e. e{p,) = {V{<py PcVitpy") ^ ^ = (^{V{<py PcVitpy") ^ 

{{vicpyp^vicpy^)^)^. 
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This same implementation idea is further generahzed to get the analogue of eq. (l37|) . For 
that s identical two-level systems are needed together with s identical and statistically in- 
dependent discrete random variables, uniformly distributed over the first k natural numbers 



^ 2w ' 



i.e. Ui ~ ({0, 1, k — 1}), i = 1,2, s, as well as an independent circular rv (p 
< (f) < 2tt). First we form the statistical average over the s independent discrete variables 

e4,{Pc) ® ... ® e^Pc) = (l^(0)VcV(0)t'^>^ ® ... ® (y(0)>,y(0)t->^ . (40) 

Then we treat the resulting s-fold tensor product of density matrices e^{pc), as (j) correlated 
matrix-valued random variables, and consider their statistical average 

® ... ® e^(p.)), = I ^iM^ (v(0)v,y(0)t-)^ ® ... ® (\/(0)v,r(0)t'^>^ . (41) 



The doubly averaged density matrix is identical to £^(pc), i.e. ^*(Pc) = 

(^^(Pc) ® ... ®£^(pc))^. 

We note finally that the range k of discrete rv ly is determined by the kind of V'^ model 
of QRW the statistical moments of which we simulate, and that similarly the number s of 
two-level atoms involved in simulation is determined by the order of quantum moment of 
walker's system we intend to simulate. Two-level atoms, and the (p, u, classical stochastic 
variables are quantum and classical resources required for this on the average stochastic 
simulation of QRW. Note also that, the main difference between quantum and stochastic 
simulations is in the way the CPTP maps of coin systems are derived: in the stochastic case 
many runs of random rotations are required so that a classical ensemble average is formed, 
while in quantal case the prescribed total unitary evolution is generated in one run and then 
an coin unconditional measurement provides the final density matrix (c.f. [34ll37j) . 

As to the experimental realization of the proposed stochastic implementation we note that 
the transformation pc V {(pY pcV {(pY" , requires a random rotation V{(f)Y = (V^/(0)[/)'^ = 
^g«</>o-3g«f o-2~)!/^ For the case of k = 2, i.e. i/ = 0,1, we need to flip randomly between the 
two rotations (l,e*'^°"^e*4'^^), of the coin system. Stochastic unitary rotations of a two-level 
atom can experimentally be achieved by randomly pulsating laser fields with appropriate 



fluctuating phases and intensities [34 1. 



As corroboration of our theoretical results about quantum simulation, we next provide 
numerical evaluation of the asymptotic probability density function of the walker. This is 
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done for the example of [/—quantization (c.f. eq.( for the models with k = 2,3), and 
the example of 5— quantization (c.f. eq. f|T9l) ). If y = h{(j)), then the cumulative probability 
function in the interval [yi , 1/2] becomes 



In the application we choose p(0, (p) = 1, corresponding to = |0)(0|. In the numerical sim- 
ulation to determine the values of the probability density function over some small interval 
/ = [a, b], we count the number of times the random variable (p, is such that yi < h{(p) < y2, 
and compute the fraction of the number of successful counts by the overall number of counts. 
The results displayed in fig. la, for the model V'^, quantized by the [/-rule, show the well 
known by now double-horn distribution (c.f. previous works cited in the introductory chap- 
ter). As to the results of fig. lb and fig. 2, they are new, and provide information about the 
asymptotic behavior of the k = 3 model quantized by [/—rule, and of the model k = 2, 
quantized according to 5— quantization rule. The results show some similarity in the form 
of the distributions, but also important differences in the displacements and the cutoffs of 
their supports. These should be detectable features in realistic simulations of a walk with 

e.ntum systems. Finally, these results should be compared and constrasted with those of 
. In that work the delayed tracing scheme with time step evolution maps ey, ^v^, ^ya, 
, was used for QRW on line, in which coin decoherence has been introduced, and modifies 
the quantum walk along the lines of e— quantization rule. These results show destruction 
of quantum features of the walk for increasing strength of decoherence. On the contrary 
results indicated in fig. 2, show that for the analogous e— quantized 5y2 model with e^a, 
n = 1,2, time step evolution maps, the quantum features are retained in the long time 
regime, albeit in a modified form. 

V. CONCLUSIONS 

Quantization rules provide a broad and systematic framework for quantizing classical 
random walks, taking into account phenomena such as losses, decoherence, noise or coher- 
ent dynamics occurring in the coin systems. Optical processes that may cause or on purpose 
induce similar phenomena in coins may therefore constitute physical probes for the study of 
novel features in quantum walks. The prospects of such quantum optical walks are further 




(42) 



yi<h{4>)<y2 
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enhanced by the possibihty of using the coin systems not only as a trigger of the walk, 
but also as a quantum simulator of its dynamics and statistics, especially of its long terms 
characteristics 35|. In this framework phenomena related to quantum-classical transitions 
in walker's evolution, and to quantum coin+walker entanglement, especially in the asymp- 
totic regime of a walk, find new theoretical and experimental possibilities for a in-depth 
investigation. To some of these topics we aim to return elsewhere. 
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programme. 
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Figure captions 
Figure la, lb 

Asymptotic probability density function, for walker's scaled position variable. It refers to 
the V'^ model in fig. la (the model in fig. lb), quantized by the U-Tule, with reshuffling 
matrix the | rotation matrix, and initial coin chosen in excited state. 

Figure 2 

Asymptotic probability density function, for walker's scaled position variable. It refers 
to the V'^ model, quantized by the e-rule, with initial coin chosen in excited state. 
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Fig. 1b U-quantized walk : simulation of asymptotic density function of the walker, k=3 




Fig. 2 Epsilon-quantized walk : simulation of asymptotic probability density function of the walker 




